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Abstract 

An extension of hydrodynamic D Angelo mode of inhomogeneous sheared plasma flow along 
the magnetic field into the short- wavelength limit, where the hydrodynamic treatment is not 
valid, has been considered. We find that D Angelo mode in this wavelength range is excited by 
inverse ion Landau damping and becomes the shear flow driven ion-kinetic mode. 

1 INTRODUCTION 

It is well established now that existence of magnetic field aligned plasma sheared flows is a fun- 
damental reality throughout the space plasmas [1, 2, 3]. The presence of velocity shear produces 
a drastic changing in the stability properties of these plasma flows, development of specific shear 
flow driven instabilities and plasma turbulence, which is responsible for efficient mechanisms for ion 
energization. A substantial amount of research has been devoted to hydrodynamic electrostatic in- 
stability of such shear flows, which was originally discussed by D'Angelo [4] and named now as a 
D'Angelo mode. This instability is of Kelvin-Helmholtz type. It was investigated experimentally and 
theoretically in different conditions specifically both in coUisionless [5] and coUisional [6] plasmas, 
in a plasmas with negative ions [7, 8] as well as in dusty plasmas [9, 10]. In homogeneous plasmas 
D'Angelo mode is a purely growing with zero frequency and growth rate 7 = kzVg {{ky/kz) S — 1)^*^^, 
where Vg — (Te/"^?) ^"^^ is the ion sound velocity, ky and kz are the components of the wave vector 
across and along the magnetic field respectively, — {l/ud) (dVo (X) /dX) is the normalized ion 
flow velocity shear, uici is the ion cyclotron frequency. Instability develops, when (ky/kz) S > 1. This 
mode has obvious difference with shear flow modified ion sound instability [11], which requires that 
(ky/kz) S < 1 and develops due to inverse electron and ion Landau damping. In Ref.[12] Mikhailenko 
et al. by using the kinetic approach, which accounted for the effects of the finite Larmor radius, elec- 
tron and ion Landau damping, have developed the comprehensive theory of the hydrodynamic and 
kinetic shear flow modified and shear flow driven low frequency ion-sound-drift instabilities in inho- 
mogeneous plasma shear flow under condition (ky/kz) S < 1- The purpose of this paper is to present 
a kinetic theory of the electrostatic low frequency instabilities of inhomogeneous plasmas under con- 
dition (ky/kz) 5" > 1, which is key to connection of the investigated instabilities with D'Angelo mode. 
In contrary to conclusions of Smith and Goeler [13], which have considered this instability on the base 
of the Vlasov equation and climed[13] that ion Landau damping leads only to modification of this 
hydrodynamical instability, we find that in short wavelength limit, where hydrodynamic treatment 
is not valid at all, the ion-kinetic effects are dominant and are responsible for the development of 



1 



new kinetic instability. Keeping in mind that this instability is characterizes by positive parameter 
(ky/kz) S, we nominate it as a kinetic D'Angelo mode. This kinetic D'Angelo mode occurs when ve- 
locity shear exceeds the critical value, which equals to one for the development of the hydrodynamic 
D'Angelo mode [4] and develops due to inverse ion Landau damping. We find, that for the velocity 
shear less or of the order of the critical value the kinetic D'Angelo mode transforms into the drift 
instability, which has a similar mechanism of excitation. 

The paper is organized as follows. In Sec. 2 the dispersion relation for D'Angelo mode 
obtained by using the kinetic approach is represented. In Sec. 3 the hydrodynamic D'Angelo mode 
for inhomogeneous plasma shortly reviewed with accounting for the effects of finite ion Larmor radius. 
Sec. 4 is devoted to analysis of the kinetic D'Angelo mode in the short wavelength limit. In Sec. 5 
the drift instability when the velocity shear less of order the critical value is considered. Also here a 
numerical analysis is presented. Conclusions are given in Sec. 6. 

2 DISPERSION RELATION FOR LOW-FREQUENCY 
INSTABILITIES DRIVEN BY SHEAR FLOW 

We consider the inhomogeneous sheared plasma flow with one ion species which moves with velocity 
Vo (x) II Boez- The general dispersion relation in the kinetic Vlasov-Poisson model was obtained in 
Ref.[14] 
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where is the Debye length, uj^a and VTa = \fTaf^a are the cyclotron frequency and 

thermal velocity respectively, An = J^n{bi)<^^^% In is the modified Bessel function, hi = 

k\pTii PTi = VTi/uJci is the thermal Larmor radius, Zan = {oj — nUca — 

{(jj — kzVoa — kyVda) /V^kzVTa-i Vda = {'^Tal^ca) (rf In tzq (x) / dx) is the diamagnctic drift velocity, and 
W {z) = exp (— ^1 + (2i/-^/7r) J exp {t^) dt^ . In Eq. (1) is assumed that the local approximation, 

kxLn I and kxLy 

^ 1, where L„ = [din no (x) /dx] , Ly — [dlnVo (x) /dx] , holds. 
We consider the low frequency instability when the condition ou — ou (k) + /c^Vq <^ oud is met, 
so the zero term in the sum over cyclotron harmonics in Eq. (1) is the dominant. We suppose that 
electrons are adiabatical with j^^eoj <S 1, and \zio\ may be greater or less than unity, whereas > 1 
for n 7^ 0. The last condition is valid, when inequality kzPn < 1 is satisfied. Using the asymptotic 
form for W - function for large argument values, W{zi) ~ exp {—z^) + (i/y/nzi) (1 + 1/22;^^), we carry 
out the summation over all cyclotron harmonics with n 7^ 0. Then dispersion relation (1) reduces to 
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the form 



k^XliS (k,cc;) = 1 + r 

+i^/^^W (zio) Aq (bi) ^Zio - XkyPnSiZiQ - XkyPTi J^' ^ - XkyPnSiAo (bi) = 0, (2) 

where r = T^/T^. In Eq. (2) the variables A = l/k^pn and Zio — {u — kzVoa) /V^kzVTa, which are 
the normahzed wavelength along the magnetic field and complex frequency respectively, are used. 
Eq. (2) is the general dispersion relation for the low frequency perturbations which takes into account 
effects of the velocity shear, thermal motion of ions, both along and across the magnetic field, as well 
as the plasma density inhomogeneity. The analytical solution to Eq. (2) in general case for arbitrary 
A and Zio values can not be found. However, it admits an approximate solution in asymptotic limits of 
long as well as short wavelengths along the magnetic field, when the effect of ion Landau damping is 
negligible or dominant, respectively. In the long wavelength limit, this equation gives hydrodynamic 
shear-fiow-driven instability, known as D'Angelo mode. In the short wavelength range, the source 
of free energy for instability is also a shear velocity, however the excitation of instability occurs via 
the inverse Landau damping and therefore this instability is nominated as the ion kinetic D'Angelo 
mode. 



3 HYDRODYNAMIC D'ANGELO MODE 

In the long wavelength limit, when inequality |2;jo| > 1 is met and ion Landau damping is negligible, 
the dispersion relation (2) reduces to the form 

24 {1 + T-Ao ih)) + V2z,oAo {hi) ^A + Ao (hi) {XkyPnSi - 1) = 0. (3) 

This equation gives the classical hydrodynamic D'Angelo mode [4], which, in addition, accounts for 
the thermal motion of ions across the magnetic field. The frequency and growth rate of this mode 
respectively are 



"(^^^y ^ = ^r kiviA,{bi)x^ {xkypr A (4) 



where 



is the drift wave frequency. It follows from Eq.(4), that the hydrodynamic D'Angelo mode occurs in 
the wavelengths range A2 < A < Ai, where 



_ 2ul,Ao ih) kyPTiSi ± ^ (2^2.^0 ih) kyPT^Sif - ujIAujI {I + t-A^ (6,)) ^0 {h) 
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The hydrodynamic D'Angelo mode exists when Ai and A2 are both real and positive. The condition 
of the wavelength positivity is satisfied when inequality kyPnSi > is met, whereas the wavelengths 
are real when the fiow velocity shear satisfies inequality 




In the limiting case k±pTi <^ 1 and t — 1, inequality (7) reduces to the condition Si > \vdi \ /vn 
obtained by D'Angelo in the fiuid treatment [4]. Under condition (7) the growth rate (4) has a 
maximum at Aq = 2/kyPTiSi. 

The condition for validation of hydrodynamic treatment, \ziQ\ > 1, at the extremities of 
interval (A2; Ai), as well as in the point of the growth rate maximum, Aq, should be verified also. At 
the long wavelength boundary, 

A = A..4..% '\^-/°'^'» . (8) 
the modulus of normalized frequency is approximately equals 

\zio\ ^ V2Si^ (9) 

\Vdi\ 

and hydrodynamic treatment is met at Si > \vdi\/ V^Vti- This inequality is the less strict than (7) 
and therefore it holds if inequality (7) is satisfied. In the point of the growth rate maximum, the 
modulus of normalized frequency is 

I I \vdi\ Ao{bi) 

and under condition (7) inequality \zio\ > 1 is not met here. This condition at A = A2 is not satisfied 
also. Hence the problem of the D'Angelo mode at the short wavelength threshold, as well as in the 
maximum of the growth rate can not be solved in the hydrodynamic treatment. 



4 KINETIC D'ANGELO MODE 

In order to investigate the D'Angelo mode in the short wave limit, when Landau damping is dominant, 
we consider Eq. (2) in general form without appealing to the asymptotic forms of W function. Find 
first the short wavelength threshold of the D'Angelo mode. The threshold values for variables A and 
Zio are determined by the balance of the velocity shear, plasma inhomogeneity as well as thermal 
motion effects in ion Landau damping and may be obtained by equating to zero the real and imaginary 
parts of Eq. (2), 

Zio - XkyPTiSiZiQ - XkyVdi I = 0, ^^^^ 
1 + r - XkyPTiSiAo (hi) = 0, 
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System of equations (11) has a solution, when inequahty kyPnSi > is met. The short-wavelength 
threshold value for the excitation of the instability, as well as the threshold value of the normalized 
complex frequency Zg, which is the real at that threshold, are 

^ 1 + r ^2 

kyPTiSiAo (bi) ' 

^ Vdi l + r 

V2S,vr.l + r-Ao{ky ^ ' 

where index s means the short- wavelength instability threshold. The approximate solution to Eq.(2) 
in the vicinity of instability threshold, we obtain by Taylor series expansion of Eq.(3) in powers of 
(A — As) with retained only zero-order and linear terms, 

Zio ~ + 4 (A,) (A - A,) . (14) 

Here z'^{Xs) — —e\/e'^ and e\ and e'^ are determined by 

k'XlA (A.) = -iV^W (z,) Ao ih) [kyPTiS.Zs + - kyPTiS.A, (6,) (15) 

exl,e', {K) = -^^/^W {z^) (1 + r - ^ (6,)) . (16) 

The expression (14) does not claim to be exact solution of dispersion relation, however it will give 
the ability to determine the condition on the velocity shear at which the D'Angelo mode in short 
wavelength range occurs. The waves frequency and the growth rate in the vicinity of the instability 
threshold can be obtained from Eq. (14) as 

A, -'o-4o(6,) / A,\ v^/,-,/;,,>5',-4o(^;)Iiiiir(.:,0 / A, 

u (k) ~ cuq- \ ' 1 ] + oJci 9 — ^ — 1 

^' A (l + r-Ao(60)V a; ''^\W{zs)\^{l + T-A,{h,))\ A 



7 - <^ci ^,^1, ,,2,. , T^TTV I 1 - — ) ' (18) 



^f2kyPTiS,A^ (bj) ReW jzs) 
\/^\W{zs)f{l + T-Ao{b,)) V A 

As follows from Eq. (18) the D'Angelo mode is a stable for the wavelengths along the magnetic field 
which are satisfied to inequality A < As, whereas for waves with A > As it becomes unstable. The 
instability is induced due to effect of inverse ion Landau damping which occurs when the free energy 
of velocity shear exceeds the absorption of energy by ion thermal motion along the magnetic field and 
plasma inhomogeneity. Formally the instability exists for arbitrary magnitudes of the velocity shear 
for wavelengths A > A^. However because 7 oc KeW (zs) = exp{—Zg), the instability growth rate 
(18) is exponentially small in the vicinity of the threshold when l^;^! > 1 and it is not exponentially 
small when l^;^! < 1. Lust inequality gives in fact the necessary condition on normalized shear, at 
which the D'Angelo mode is unstable. Accounting for Eq.(13) yields 

^^^^^°^l + r-Io(&OV^l- ^'^^ 
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In the kyPTi <^ 1 limit and for r = 1 inequality (19) reduces to Si > \/2\vdi \ /vti which coincides 
with that obtained by D'Angelo in the fluid treatment [4]. The magnitude of the growth rate is 
affected also by factor kyPTiA^ / {1 + t — Aq which gives the maximum of the growth rate at 
kyPTi - 0.5. 

This kinetic instability reveals the main properties of the hydrodynamic D'Angelo mode, such 
as an increase in growth rate with increasing of velocity shear and it decreases with increasing of 
plasma inhomogeneity. The latter property, however, is manifested only in the quadratic term of 
the expansion (14), which is omitted here because of its lengthy form. It is interesting to note, that 
the boundary wavelength (12) for the kinetic mode is independent on the plasma inhomogeneity, in 
distinction from hydrodynamic mode. Thus, the D'Angelo mode in the vicinity of the short wave 
threshold occurs due to inverse ion Landau damping which is caused by the parallel velocity shear 
and it is the ion-kinetic D'Angelo mode. The necessary condition on velocity shear value is expressed 
by Eq. (19). The most unstable waves are those with transverse wave numbers kyPTi — 0.5. 

The growth rate of the ion-kinetic D'Angelo mode (18) increases with an increase of wave- 
length until condition \ziQ\ < 1 holds. With subsequent increases of the wavelength, so that inequality 
l^iol > 1 is met, the role of kinetic effects on the development of the D'Angelo mode reduces, whereas 
the hydrodynamic mechanism becomes dominant and ion-kinetic D'Angelo mode transformes grad- 
ually into the classical hydrodynamic D'Angelo mode. In the transition wavelength region, where 
\zio\ ^ 1 and both the kinetic and hydrodynamic effects are significant, an analytic solution of Eq. (2) 
not obtained. The corresponding range of the wavelengths along the magnetic field may be studied 
only numerically. 



5 ION-KINETIC SHEAR-FLOW-DRIVEN DRIFT INSTA- 
BILITY 

When the velocity shear satisfies inequality Si < Sio, so that Zs > 1, the D'Angelo mode does not 
exist, whereas the kinetic growth rate, though small compared to D'Angelo mode, but is not equal 
to zero. In the asymptotic limit Zio > Zg > 1, the dispersion relation (2) for inhomogeneous plasmas 
has a form 

24 {1 + T-Ao ih)) + V2zioAo {hi) + Ao (k) {Xkyp^Si - 1) 

^-li^fnAQ (bi) 4 (^Zio - XkyPnSiZiQ - XkyPn J exp (-4) = 0, (20) 

which is almost a same as Eq. (3), but we retain in the (20) the ion Landau damping term. The 
frequency of instability, 

LJ (k) ^ LJo + kzVTi—Si, (21) 
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Figure 1: The growth rate of the D'Angelo mode vs the normahzed wavelength: (1) Si — 0.3, (2) 
Si = 0.2, (3) Si = 0.1. Here \vdi\/vTi = 0.1, kyPn = 0.5, r = 1. 

corresponds to the drift waves modified by velocity shear. Therefore this waves may be identified as 
a drift waves. The growth rate of the drift instability. 
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(,^,„SA - ^) , ^^''l' exp (-'^) (22) 




is positive for wavelengths A > Ag, so, the drift instability, as well as the ion kinetic d'Angelo mode, is 
excited due to inverse ion Landau damping. The maximum of the growth rate occurs approximately 
at wavelength 

\ / / OC2\ 

(23) 

With increasing wavelengths at A > Aq, the growth rate decreases exponentially. 

In order to verify the results of analytical calculations, we solved numerically the dispersion 
relation (1) in the frequency range uu <^ ouci- Figure 1 shows the growth rate of instability versus the 
normalized wavelength for the different magnitudes of normalized velocity shear. Si — 0.3, Si = 0.2, 
and Si = 0.1 (curves 1, 2 and 3), which appropriate to ionospheric parameters [1]. Here we take also 
\vdi\/vTi = 0.1, kyPTi = 0.5, r = 1. The condition for velocity shear (19) is satisfied only for the first 
two values of shear, so that only curves 1 and 2 present the D'Angelo mode, whereas curve 3 presents 
the drift instability. For these conditions, the value of l^;^! at the short wavelength boundary is equal 
to 0.4, 0.6 and 1.2, respectively. The short-wavelength boundary A^ calculated from Eq. (13) is equal 
approximately to 17, 25 and 51, respectively, which agrees with Figure data. The long- wavelength 
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boundary Ai of the hydrodynamic D'Angelo mode for Si — 0.3 and Si — 0.2 was calculated from Eq. 
(9) and is equal to 360 and 235, respectively. The growth rate for curves 1 and 2 has a maximum 
at wavelengths corresponding to l^jol — 1 i- e. in the transient kinetic-hydrodynamic range of the 
longitudinal wavelengths, where the kinetic and hydrodynamic effects are significant both. For shear 
value Si — 0.1 (curve 3), the condition of hydrodynamic approach (8) does not satisfied, and thus 
D'Angelo mode does not develop. So this mode is the kinetic shear-flow driven drift instabihty. The 
point of the maximal growth rate for this instability, calculated from Eq. (23), is approximately 
equal A 65, which is shght less than the value from the Figure. 

6 CONCLUSIONS 

An extension of hydrodynamic D'Angelo mode of inhomogcncous sheared plasma flow along the mag- 
netic field into the short-wavelength limit, where the hydrodynamic treatment is not valid, has been 
considered. We find that D'Angelo mode in this wavelength range is excited by inverse ion Landau 
damping and becomes the shear fiow driven ion-kinetic mode. The short-wavelength boundary for 
the wave number component along the magnetic field A^, (12) of instability development is determined 
by the balance of shear fiow, plasma inhomogeneity and thermal motion effects in ion Landau damp- 
ing. The kinetic D'Angelo mode is induced for the wavelengths A > A^ and velocity shear S > SiQ, 
where Sio = {\vdi\ /\/2vTi) {I + r) / {1 + t — Aq (bi)) is the critical value. This condition for velocity 
shear is the same as for the ordinary D'Angelo mode [4], so that kinetic and hydrodynamic D'Angelo 
modes exist simultaneously both in different wavelength regions. The kinetic D'Angelo mode reveals 
the main properties of the hydrodynamic D'Angelo mode, such as an increase in growth rate with 
increasing of velocity shear; and it decreases with increasing of plasma inhomogeneity. However, 
the short wavelength boundary (11) for the kinetic mode appears to be independent on the plasma 
inhomogeneity unlike the long wavelength boundary (8) of hydrodynamic mode. The maximum of 
the growth rate of D'Angelo mode occurs in transient wavelength range, where both the kinetic and 
hydrodynamic effects are significant. For velocity shear value S < -S'io and wavelengths A < As the 
D'Angelo mode does not exist, however, the ion-kinetic drift instability at A > As is developed due to 
velocity shear. This instability occurs as the kinetic D'Angelo mode, when the free energy of velocity 
shear exceeds the absorption of energy by ion thermal motion along the magnetic field and plasma 
inhomogeneity. Formally, this instability exists for arbitrary magnitudes of the velocity shear for 
wavelengths A > As, however its growth rate (19) is exponentially small, when \zs\ ^ 1. 
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